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Abstract I begin with a proposed global phase diagram of the cuprate supercon- 
ductors as a function of carrier concentration, magnetic field, and temperature, and 
highlight its connection to numerous recent experiments. The phase diagram is then 
used as a point of departure for a pedagogical review of various quantum phases 
and phase transitions of insulators, superconductors, and metals. The bond opera- 
tor method is used to describe the transition of dimerized antiferromagnetic insula- 
tors between magnetically ordered states and spin-gap states. The Schwinger boson 
method is applied to frustrated square lattice antiferromagnets: phase diagrams con- 
taining collinear and spirally ordered magnetic states, Z2 spin liquids, and valence 
bond solids are presented, and described by an effective gauge theory of spinons. 
Insights from these theories of insulators are then applied to a variety of symme- 
try breaking transitions in c?-wave superconductors. The latter systems also contain 
fermionic quasiparticles with a massless Dirac spectrum, and their influence on the 
order parameter fluctuations and quantum criticality is carefully discussed. I con- 
clude with an introduction to strong coupling problems associated with symmetry 
breaking transitions in two-dimensional metals, where the order parameter fluctua- 
tions couple to a gapless line of fermionic excitations along the Fermi surface. 



1 Introduction 

The cuprate superconductors have stimulated a great deal of innovative theoretical 
work on correlated electron systems. On the experimental side, new experimental 
techniques continue to be discovered and refined, leading to striking advances over 
twenty years after the original discovery of high temperature superconductivity |l]|. 
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Fig. 1 Proposed global phase diagram for the hole- and electron-doped cuprates [3 -51. The axes 
are carrier concentration (x), temperature (T), and a magnetic field (H) applied perpendicular to 
the square lattice. The regions with the SC label have d-wave superconductivity. The strange metal 
and the "pseudogap" regime are separated by the temperature T* . Dashed lines indicate crossovers. 
After accounting for the valence bond solid (VBS) or Ising nematic orders that can appear in the 
regime Xg < x < Xm, the dashed T* line and the dotted line connecting Xm to the point M 
become true phase transitions. There can also be fractionalized phases in the region Xs < x < Xm 
, as discussed recently in Refs. [^jTj. 

In the past few years, a number of experiments, and most especially the discovery 
of quantum oscillations in the underdoped regime \2j, have shed remarkable new 
light on the origins of cuprate superconductivity. I believe these new experiments 
point to a synthesis of various theoretical ideas, and that a global theory of the rich 
cuprate phenomenology may finally be emerging. The ingredients for this synthesis 
were described in Refs. |[3]-Q, and are encapsulated in the phase diagram shown in 
Fig.[T] Here I will only highlight a few important features of this phase diagram, and 
use those as motivations for the theoretical models described in these lectures. The 
reader is referred to the earlier papers ||4]|5) for a fuU discussion of the experimental 
support for these ideas. Throughout the lectures, I will refer back to Fig.[T]and point 
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out the relevance of various field theories to different aspects of this rich phase 
diagram. 

It is simplest to examine the structure of Fig.[T]beginning from the regime of large 
doping. There, ample evidence has established that the ground state is a conventional 
Fermi liquid, with a single "large" Fermi surface enclosing the area demanded by 
the Luttinger theorem. Because of the underlying band structure, this large Fermi 
surface is hole-like (for both electron and hole-doping), and so encloses an area 
1+x for hole density x, and an area 1— p for doped electron density p. The central 
quantum phase transition (QPT) in Fig.[T]is the onset of spin density wave (SDW) 
order in this large Fermi surface metal at carrier concentration x ~ Xm, shown in 
Fig. [T] near the region where Tc is largest (the subscript m refers to the fact that 
the transition takes place in a metal); we will describe this transition in more detail 
in Section]?] Because of the onset of superconductivity, the QPT at x = x,n is 
revealed only at magnetic fields strong enough to suppress superconductivity, i.e., at 
H > Hc2- For X < Xm, we then have a Fermi liquid metal with SDW order. Close 
to the transition, when the SDW order is weak, the large Fermi surface is generically 
broken up by the SDW order into "small" electron and hole pockets, each enclosing 
an area of order x (see Fig. 17 later in the text). Note that electron pockets are present 
both for hole and electron doping: such electron pockets in the hole-doped cuprates 
were first discussed in Ref. |[8|. There is now convincing experimental evidence for 
the small Fermi pockets in the hole underdoped cuprates, including accumulating 
evidence for electron pockets ||9][T0). The QPT between the small and large Fermi 
surface metals is believed to be at Xm ~ 0.24 in the hole-doped cuprates |jTT 12 1, 
and atpm ~ 0.165 in the electron-doped cuprates | T3p4 |. One of the central claims 
of Fig.[T]is that it is the QPT at a; = Xm which controls the non-Fermi liquid "strange 
metal" behavior in the normal state above the superconductivity Tc- We leave open 
the possibility |15 | that there is an extended non-Fermi liquid phase for a range of 
densities with x > Xm- this is not shown in Fig.[T] and will be discussed here only 
in passing. 

The onset of superconductivity near the SDW ordering transition of a metal has 
been considered in numerous previous works p6p7 ) . These early works begin with 
the large Fermi surface found for x > Xm, and consider pairing induced by ex- 
change of SDW fluctuations; for the cuprate Fermi surface geometry, they find an 
attractive interaction in the d-wave channel, leading to d-wave superconductivity. 
Because the pairing strength is proportional to the SDW fluctuations, and the latter 
increase as a; \ Xm, we expect Tc to increase as x is decreased for x > Xm, as 
is shown in Fig. [T] Thus for x > Xm, stronger SDW fluctuations imply stronger 
superconductivity, and the orders effectively attract each other. 

It was argued in Ref. |3J that the situation becomes qualitatively different for x < 
Xm- This becomes clear from an examination of Fig. [T] as a function of decreasing 
T for X < Xm- It is proposed |3, 18| in the figure that the Fermi surface already 
breaks apart locally into the small pocket Fermi surfaces for T < T*. So the onset 
of superconductivity at Tc involves the pairing of these small Fermi surfaces, unlike 
the large Fermi surface pairing considered above for x > Xm- For x < Xm, an 
increase in local SDW ordering is not conducive to stronger superconductivity: the 
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SDW order 'eats up' the Fermi surface, leaving less room for the Cooper pairing 
instability on the Fermi surface. Thus in this regime we find a competition between 
SDW ordering and superconductivity for 'real estate' on the Fermi surface |[3 19 1. 
As we expect the SDW ordering to increase as x is decreased for x < Xm, we should 
have a decrease in Tc with decreasing x, as is indicated in Fig.[T] 

We are now ready to describe the second important feature of Fig.[T] The comple- 
ment of the suppression of superconductivity by SDW ordering is the suppression 
of SDW ordering by superconductivity. The competition between superconductivity 
and SDW order moves |3| the actual SDW onset at H = and T = to a lower 
carrier concentration x = Xg (ot p — ps for electron doping). The QPT at a; = Xg 
controls the criticality of spin fluctuations within the superconducting phase (and 
hence the subscript s), while that x — Xm continues to be important for T > Tc 
(as is indicated in Fig. [TJ. There is now a line of SDW-onset transitions within 
the superconducting phase pUj , connecting the point Xs to the point M, for which 
there is substantial experimental evidence |21fj26|. The magnitude of the shift from 
Xm to Xs depends a great deal upon the particular cuprate: it is largest in the ma- 
terials with the strongest superconductivity and the highest Tc. In the hole-doped 
YBCO series we estimate Xg ~ 0.085 \26] and in the hole-doped LSCO series 
we have Xg ~ 0.14 |23| (recall our earlier estimate Xm ~ 0.24 in the hole-doped 
cuprates 1 11 12)), while in the electron-doped cuprate Nd2_xCea;Cu04, we have 
Ps = 0.145 f21] (recall p„i ~ 0.165 in the electron-doped cuprates P131I14|). 

With the shift in SDW ordering from Xm to Xg, the need for the crossover line la- 
beled Tsdw in Fig.[T|becomes evident. This is the temperature at which the electrons 
finally realize that they are to the 'disordered' of the actual SDW ordering transition 
at a; = Xs, rather than to the 'ordered' side of the transition at a; = x„i- Thus, for 
T < Tsdw, the large Fermi surface re-emerges at the lowest energy scales, and SDW 
order is never established. This leaves us with an interesting superconducting state 
at T = 0, where the proximity to the Mott insulator can play an important role. 
Other orders linked to the antiferromagnetism of the Mott insulator can appear here, 
such as valence bond solid (VBS) and Ising-nematic order |28|, or even topologi- 
cally ordered phases ||6||7]; experimental evidence for such orders has appeared in 
a number of recent experiments p2l[29||3T[ , and we will study these orders in the 
sections below. 

The shift in the SDW ordering from x,n to Xs has recently emerged as a generic 
property of quasi-two-dimensional correlated electron superconductors, and is not 
special to the cuprates. Knebel et al. |32| have presented a phase diagram for 
CeRhlns as a function of temperature, field, and pressure (which replaces carrier 
concentration) which is shown in Fig.|2] Notice the very similar structure to Fig.[T] 
the critical pressure for the onset of antiferromagnetism shifts from the metal to the 
superconductor, so that the range of antiferromagnetism is smaller in the supercon- 
ducting state. In the pnictides, the striking observations by the Ames group |33 ■ 35) 
on Ba[Fei_2;Co2;]2As2 show a 'back-bending' in the SDW onset temperature upon 
entering the superconducting phase: see Fig. [2] This is similar to the back-bending 
of the line Tgdw from T* in Fig.[T[ and so can also be linked to the shift in the SDW 
onset transition between the metal and the superconductor. 
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Fig. 2 Phase diagrams for CeRhlns from Ref. (32| and for Ba[Fei_i,Coj:]2As2 from Refs. j33j- 
[35 1. For CeRhlns, the shift from pc to p* is similar to the shift from Xm to Xs in Fig.[T] this 
shift is significantly larger in the cuprates (and especially in YBCO) because the superconductiv- 
ity is stronger. In the ferropnictide Ba[Fei_a:Coi,]2As2, the back-bending of the SDW ordering 
transition in the superconducting phase is similar to that of Tg^^ in Fig.[T| 



It is clear from Fig.|2]that the shift in the SDW order between the metal and the 
superconductor is relatively small in the non-cuprate materials, and may be over- 
looked in an initial study without serious consequences. Similar comments apply 
to the electron-doped cuprates. However, the shift is quite large in the hole-doped 
cuprates: this can initially suggest that the cuprates are a different class of materials, 
with SDW ordering playing a minor role in the physics of the superconductivity. 
One of the main claims of Fig. [T] is that after accounting for the larger shift in 
the SDW transition, all of the cuprates fall into a much wider class of correlated 
electron superconductors for which the SDW ordering transition in the metal is the 
central QPT controlling the entire phase diagram (see also the recent discussion by 
Scalapino |36|). 

Our discussion will be divided into 3 sections, dealing with the nature of quantum 
fluctuations near SDW ordering in insulators, d-wave superconductors, and metals 
respectively. These cases are classified according to the increasing density of states 
for single-electron excitations. We will begin in Section |2] by considering a vari- 
ety of Mott insulators, and describe their phase diagrams. The results will apply 
directly to experiments on insulators not part of the cuprate family. However, we 
will also gain insights, which will eventually be applied to various aspects of Fig.[T] 
for the cuprates. Then we will turn in Section|3]to d-wave superconductors, which 
have a Dirac spectrum of single-electron excitations as described in Section |3.1| 
Their influence on the SDW ordering transition at a; = a;^ will be described using 
field-theoretical methods in Section 3.2 Section [33] will describe the Ising-nematic 
ordering at or near x — Xg indicated in Fig.[T] Finally, in Section]?] we will turn 
to metals, which have a Fermi surface of low-energy single-particle excitations. We 
will summarize the current status of QPTs of metals: in two dimensions most QPTs 
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lead to strong coupling problems which have not been conquered. It is clear from 
Fig. [TJthat such QPTs are of vital importance to the physics of metallic states near 



Significant portions of the discussion in the sections below have been adapted 
from other review articles by the author 1 38 -401. 



2 Insulators 

The insulating state of the cuprates at x = is a 5 = 1/2 square lattice antifer- 
romagnet, which is known to have long-range Neel order. We now wish to study 
various routes by which quantum fluctuations may destroy the Neel order. In this 
section, we will do this by working with undoped insulators in which we modify 
the exchange interactions. These do not precisely map to any of the transitions in 
the phase diagram in Fig.[T] but we will see in the subsequent sections that closely 
related theories do play an important role. 

The following subsections will discuss two distinct routes to the destruction of 



Neel order in two-dimensional antiferromagnet. In Section 2.1 we describe coupled 
dimer antiferromagnets, in which the lattice has a natural dimerized structure, with 
2 5* = 1/2 spins per unit cell which can pair with each other. These are directly 
relevant to experiments on materials like TlCuCla. We will show that these anti- 
ferromagnets can be efficiently described by a bond-operator method. Then in sec- 
tion 2.2 we will consider the far more complicated and subtle case where the lattice 
has full square lattice symmetry with only a single S" = 1 /2 spin per unit cell, and 
the Neel order is disrupted by frustrating exchange interactions. We will explore the 
phase diagram of such antiferromagnets using the Schwinger boson method. These 
results have direct application to experimental and numerical studies of a variety 
of two-dimensional Mott insulators on the square, triangular, and kagome lattices; 
such applications have been comprehensively reviewed in another recent article by 
the author [37^ , and so will not be repeated here. 



2.1 Coupled dimer antiferromagnets: bond operators 

We consider the "coupled dimer" Hamiltonian pT) 

Hd = J J2 S,-Sj+gJ S,-Sj, (1) 

where Sj are spin- 1/2 operators on the sites of the coupled-ladder lattice shown 
in Fig. [3] with the A links forming decoupled dimers while the B links couple the 
dimers as shown. The ground state of Hd depends only on the dimensionless cou- 
pling g, and we will describe the low temperature (T) properties as a function of g. 
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Fig. 3 The coupled dimer antiferromagnet. Spins (5 = 1/2) are placed on the sites, the A links 
are shown as full lines, and the B links as dashed lines. 



We will restrict our attention to J > and < 5 < 1. A three-dimensional model 
with the same structure as Hd describes the insulator TlCuCla p2f|44| . 

Note that exactly at g = 1, Hd is identical to the square lattice antiferromagnet, 
and this is the only point at which the Hamiltonian has only one spin per unit cell. 
At all other values of g, Hd has a pair of S* = 1/2 spins in each unit cell of the 
lattice. 



2.1.1 Phases and their excitations 

Let us first consider the case where g is close to 1. Exactly at g = 1, Hd is identical 
to the square lattice Heisenberg antiferromagnet, and this is known to have long- 
range, magnetic Neel order in its ground state, i.e., the spin-rotation symmetry is 
broken and the spins have a non-zero, staggered, expectation value in the ground 
state with 

(S,) = 7j,Non, (2) 

where n is some fixed unit vector in spin space, rjj is ±1 on the two sublattices, and 
A'o is the Neel order parameter This long-range order is expected to be preserved for 
a finite range of g close to 1 . The low-lying excitations above the ground state consist 
of slow spatial deformations in the orientation n: these are the familiar spin waves, 
and they can carry arbitrarily low energy, i.e., the phase is 'gapless' . The spectrum of 
the spin waves can be obtained from a text-book analysis of small fluctuations about 
the ordered Neel state using the Holstein-Primakoff method |45|: such an analysis 
yields two polarizations of spin waves at each wavevector k = {kx, ky) (measured 
from the antiferromagnetic ordering wavevector), and they have excitation energy 
= (c^^x + "^y^y)^^^' ^''■^ Cj.,Cj, the spin-wave velocities in the two spatial 
directions. 

Let us turn now to the vicinity of g = 0. Exactly at g = 0, Hd is the Hamilto- 
nian of a set of decoupled dimers, with the simple exact ground state wavefunction 
shown in Fig.|4j the spins in each dimer pair into valence bond singlets, leading to 
a paramagnetic state which preserves spin rotation invariance and all lattice sym- 
metries. Excitations are now formed by breaking a valence bond, which leads to a 
three-fold degenerate state with total spin 5 = 1, as shown in Fig. |5^. At g = 0, 
this broken bond is localized, but at finite g it can hop from site-to-site, leading to 
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Fig. 4 Schematic of the quantum paramagnet ground state for small g. The ovals represent singlet 
valence bond pairs. 




■ - 4 



(a) 



(b) 



Fig. 5 (a) Cartoon picture of the bosonic S = 1 excitation of the paramagnet. (b) Fission of the 
S = 1 excitation into two S = 1/2 spinous. The spinous are connected by a "string" of valence 
bonds (denoted by dashed ovals) which lie on weaker bonds; this string costs a finite energy per 
unit length and leads to the confinement of spinous. 



a triplet quasiparticle excitation. Note that this quasiparticle is not a spin-wave (or 
equivalently, a 'magnon') but is more properly referred to as a spin 1 exciton or 
a triplon. We parameterize its energy at small wavevectors k (measured from the 
minimum of the spectrum in the Brillouin zone) by 

ek = A+ y \ (3) 

where A is the spin gap, and c^, Cy are velocities; we will provide an explicit deriva- 
tion of (|3]l in Section p.l.2| Fig.|5]also presents a simple argument which shows that 
the S = 1 exciton cannot fission into two 5 = 1/2 'spinons'. 

The very distinct symmetry signatures of the ground states and excitations be- 
tween g ~ 1 and .g ~ make it clear that the two limits cannot be continuously 
connected. It is known that there is an intermediate second-order phase transition 
at [41, 46] g = Oc = 0.52337(3) between these states as shown in Fig.|6] Both the 
spin gap A and the Neel order parameter A'o vanish continuously as gc is approached 
from either side. 
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Fig. 6 Ground states of Hj^ as a function of g. The quantum critical point is at |46| gc = 
0.52337(3). The compound TlCuCla undergoes a similar quantum phase transition under applied 
pressure |42||44| . 



2.1.2 Bond operators and quantum field theory 

In this section we will develop a continuum description of the low energy excitations 
in the vicinity of the critical point postulated above. There are a number of ways to 
obtain the same final theory: here we will use the method of bond operators r47"48|, 
which has the advantage of making the connection to the lattice degrees of freedom 
most direct. We rewrite the Hamiltonian using bosonic operators which reside on 
the centers of the A links so that it is explicitly diagonal at g — 0. There are 4 states 
on each A link (|tt)' Iti)^ lit)' ^nd and we associate these with the canonical 
singlet boson s and the canonical triplet bosons ta (a = x, y, z) so that 

\s)^s^\Q)^^{\n)-\m ; iu^4io> = ^(in>-u;» ; 
it,)^4io) = -^(in> + u;)) ; |i.)^4lo) = ^(in> + ut». (4) 

Here |0) is some reference vacuum state which does not correspond to a physical 
state of the spin system. The physical states always have a single bond boson and so 
satisfy the constraint 

S^S + tlta = 1 . (5) 

By considering the various matrix elements (s|5'i|ia), {s\S2\ta), . . ., of the spin 
operators Si 2 on the ends of the link, it follows that the action of Si and ^2 on the 
singlet and triplet states is equivalent to the operator identities 




10 



Subir Sachdev 



S2a = ^ (^-sUa - tls - ieabctlt^ ■ (6) 

where a,h,c take the values x,y,z, repeated indices are summed over and e is the to- 
tally antisymmetric tensor Inserting (|6]) into Q, and using (|5]l, we find the following 
Hamiltonian for the bond bosons; 



Hd — Ho + Hi 

■4".- ■ 4 



3 1 

Ho^ jJ2 (~j4-''i + jAjta 



'meA 



a{£, m) {tljynaS^m^t + 4a4ia«mSf + H.C.^ + b{i, m) 



tlljtecSm, + H.C.^ + c{i,m) (tl^tl^^t^nbttb — ^la^mb^rnattb 



(7) 



where £, m label links in A, and a, 6, c are numbers associated with the lattice cou- 
plings which we will not write out explicitly. Note that Hi = at g = 0, and so 
the spectrum of the paramagnetic state is fully and exactly determined. The main 
advantage of the present approach is that application of the standard methods of 
many body theory to (j7]i, while imposing the constraint (jSj, gives a very satisfactory 
description of the phases with g ^ {), including across the transition to the Neel 
state. In particular, an important feature of the bond operator approach is that the 
simplest mean field theory already yields ground states and excitations with the cor- 
rect quantum numbers; so a strong fluctuation analysis is not needed to capture the 
proper physics. 

A complete numerical analysis of the properties of (j7]l in a self-consistent 
Hartree-Fock treatment of the four boson terms in Hi has been presented in 
Ref. pT) . In all phases the s boson is well condensed at zero momentum, and the 
important physics can be easily understood by examining the structure of the low- 
energy action for the ta bosons. For the particular Hamiltonian ([T]l, the spectrum of 
the ta bosons has a minimum at the momentum (0, tt), and for large enough g the ta 
condense at this wavevector: the representation (j6]l shows that this condensed state 
is the expected Neel state, with the magnetic moment oscillating as in Q. The con- 
densation transition of the ta is therefore the quantum phase transition between the 
paramagnetic and Neel phases of the coupled dimer antiferromagnet. In the vicinity 
of this critical point, we can expand the ta bose field in gradients away from the 
(0, tt) wavevector: so we parameterize 

iAa(T)=ta(r,,T)e'("'^)"-^ (8) 



where r is imaginary time, r = (x, y) is a continuum spatial coordinate, and expand 
the effective action in spatial gradients. In this manner we obtain 
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Sf = / d^rdr 



4^ + Ctlta - ^ {tJa + H.C.) + Ki^ld.Jal'' + Kiy\dyta\ 
1 



+ 2 iK2AdJaf + K2y{dytaf + H.C.) + 



(9) 



Here C, D, Ki,2x.y are constants that are determined by the solution of the self- 
consistent equations, and the ellipses represent terms quartic in the ta- The action 
St can be easily diagonahzed, and we obtain a 5* = 1 quasiparticle excitation with 
the spectrum 



(D + K^xk 



1/2 



(10) 



This is, of course, the triplon (or spin exciton) excitation of the paramagnetic phase 
postulated e arlier in (|3] l; the latter result is obtained by expanding ( 10 1 in momenta, 
with A — ^JC^ — D"^. This value of A shows that the ground state is paramagnetic 
as long as C > D, and the quantum critical point to the Neel state is at C = D. 

The critical point and the Neel state are more conveniently described by an al- 
ternative formulation of St (although an analysis using bond operators directly is 
also possible |49J). It is useful to decompose the complex field ta into its real and 
imaginary parts as follows 

ta= Z{Lpa + i'^a), (H) 

where Z is a normalization chosen below. From ([H]) and the connection to the lattice 
spin operators, it is not difficult to show that the vector tp^ is proportional to the 
Neel order parameter n in Eq. (j2]). Insertion of ( 11 1 into (j9]l shows that the field iTa 
has a quadratic term ~ (C + D)^^, and so the coefficient of tt^ remains large even 
as the spin gap A becomes small. Consequently, we can safely integrate tTq out, and 
the resulting action for the Neel order parameter cpa takes the form 



S^ 



d rdr 



(12) 

Here we have chosen Z to fix the coefficient of the temporal gradient term, and 
The action S^p gives a simple picture of excitations across the quantum critical 



point, which can be quantitatively compared to neutron scattering experiments |^ 
on TICUCI3. In the paramagnetic phase (s > 0), a triplet of gapped excitations is 
observed, corresponding to the three normal modes of ipa oscillating about (pa ~ 0; 
as expected, this triplet gap vanishes upon approaching the quantum critical point. In 
a mean field analysis, the field theory in Eq. ( 12 1 has a triplet gap of -/s (mean field 
theory is applicable to TlCuCls because this antiferromagnet is three dimensional). 
In the Neel phase, the neutron scattering detects 2 gapless spin waves, and one 
gapped longitudinal mode pO] (the gap to this longitudinal mode vanishes at the 
quantum critical point), as is expected from fluctuations in the inverted 'Mexican 
hat' potential of S^p for s < 0. The longitudinal mode has a mean-field energy 
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Fig. 7 Energies of the gapped collective modes across the pressure (p) tuned quantum phase transi- 
tion in TlCuCls observed by Ruegg et al. 1 44 1 . We test the description by the action S^p in Eq. ^12^ 
with s oc (pc — p) by comparing \/2 times the energy gap for p < pc with the energy of the 
longitudinal mode for p > pc. The lines are the fits to a -^Ip — pc| dependence, testing the 1/2 
exponent. 



gap of ■y/2|s|. These mean field predictions for the energy of the gapped modes 
on the two sides of the transition are tested in Fig. |7 the observations are in good 
agreement with the 1/2 exponent and the predicted 1 37 sTj \/2 ratio, providing a 
non-trival experimental test of the field theory. 

We close this subsection by noting that all of the above results have a direct gen- 
erahzation to other lattices. One important difference that emerges in such calcula- 
tions on some frustrated lattices |52| is worth noting explicitly here: the minimum 
of the ta spectrum need not be at special wavevector like (0, tt), but can be at a more 
generic wavevector K such that K and —K are not separated by a reciprocal lattice 
vector. A simple example which we consider here is an extension of ([T]l in which 
there are additional exchange interactions along all diagonal bonds oriented 'north- 
east' (so that the lattice has the connectivity of a triangular lattice). In such cases, 
the structure of the low energy action is different, as is the nature of the magnetically 
ordered state. The parameterization dSll must be replaced by 



(13) 



where 2a are independent complex fields. Proceeding as above, we find that the 



low-energy effective action ( 12 1 is replaced by 



(P'rdr 
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\(1> 



(14) 



where now <Pa is a complex field such that {'Pa) ^ {tia) ^ (4a)- Notice that 
there is now a second quartic term with coefficient v.lf v > 0, configurations with 
#Q = are preferred: in such configurations <Pa — nia + in2a, where ni.2a are two 
equal-length orthogonal vectors. Then from ( 13 1 and (|6]l it is easy to see that the 



physical spins possess spiral order in the magnetically ordered state in which is 
condensed. For the case t; < 0, the optimum configuration has (Pa — riae^^ where 
Ua is a real vector: this leads to a magnetically ordered state with spins polarized 
collinearly in a spin density wave at the wavevector K. The critical properties of 
the model in Eq. 



14 1 have been described in Ref. ||53J. 



2.2 Frustrated square lattice antiferromagnets: Schwinger bosons 

As discussed at the beginning of Section [2j the more important and complex cases 
of quantum antiferromagnets are associated with those that have a single S ^ 1/2 
spin per unit cell. Such models are more likely to have phases in which the exotic 
spinon excitations of Fig. |5]are deconfined, i.e., their ground states possess neutral 
S = 1/2 excitations and 'topological' order. We will meet the earliest established 
examples ]54||55| of such phases below. 

We are interested in Hamiltonians of the form 

H-^JyS, -5, (15) 

where we consider the general case of Si being spin S quantum spin operators on 
the sites, i, of a 2-dimensional lattice. The are short-ranged antiferromagnetic 
exchange interactions. We will mainly consider here the so-called square lattice Ji- 
J2-J3 model, which has first, second, and third neighbor interactions (see Fig. IS]). 
Similar results have also been obtained on the triangular and kagome lattices 1561 

The main direct applications of the results here are to experiments on a variety 
of two-dimensional Mott insulators on the square, triangular, and kagome lattices. 
As noted earlier, we direct the reader to Ref. |37 | for a discussion of these ex- 
periments. There have also been extensive numerical studies, also reviewed in the 



previous article |37|, which are in good accord with the phase diagrams presented 
below. Applications to the cuprates, and to Fig.[T[ will be discussed in the following 
sections. 

A careful examination of the non-magnetic 'spin-liquid' phases requires an ap- 
proach which is designed explicitly to be valid in a region well separated from Neel 
long range order, and preserves SU(2) symmetry at all stages. It should also be de- 
signed to naturally allow for neutral S = 1/2 excitations. To this end, we introduce 
the Schwinger boson description |58 |, in terms of elementary S ~ 1/2 bosons. For 
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Fig. 8 The J1-J2-J3 antiferromagnet. Spin S spins are placed on each site of the square lattice, 
and they are coupled to all first, second, and third neighbors as shown. The Hamiltonian has the 
full space group symmetry of the square lattice, and there is only one spin per unit cell. 



the group SU{2) the complete set of (25 + 1) states on site i are represented as 
follows 

\S,m) = ^ ^(6U^+"(6M^-"'|0), (16) 

where m ~ ~S, ... 5* is the z component of the spin (2m is an integer). We have 
introduced two flavors of bosons on each site, created by the canonical operator b]^, 
with a =ti U ™d |0) is the vacuum with no bosons. The total number of bosons, rib 
is the same for all the states; therefore 

bib? = ub (17) 

with Tib = 25 (we will henceforth assume an implied summation over repeated 
upper and lower indices). It is not difficult to see that the above representation of the 
states is completely equivalent to the following operator identity between the spin 
and boson operators 

2 

where a — x,y, z and the cr" are the usual 2x2 Pauli matrices. The spin-states on 
two sites i, j can combine to form a singlet in a unique manner - the wavefunction 
of the singlet state is particularly simple in the boson formulation: 



S:a = ^bla^%b^ , (18) 



{^''^blbl) |0) . (19) 

Finally we note that, using the constraint ( [T7] i, the following Fierz-type identity can 
be established 

{^"^blb]^) {e^sb]b'^) = -2S, ■ + nl/2 + S,,nb (20) 
where e is the totally antisymmetric 2x2 tensor 



1 

-1 



(21) 
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This implies that % can be rewritten in the form (apart from an additive constant) 

^ = - ^ E (^"^^l^l/j) (^7^^7&|) ■ (22) 

This form makes it clear that Ti. counts the number of singlet bonds. 

We have so far defined a one-parameter {rib) family of models Ti. for a fixed 
realization of the Jij. Increasing ni, makes the system more classical and a large 
rih expansion is therefore not suitable for studying the quantum-disordered phase. 
For this reason we introduce a second parameter - the flavor index a on the bosons 
is allowed to run from 1 . . . 2N with N an arbitrary integer This therefore allows 
the bosons to transform under SU{2N) rotations. However the SU {2N) symmetry 
turns out to be too large. We want to impose the additional restriction that the spins 
on a pair of sites be able to combine to form a singlet state, thus generalizing the 
valence-bond structure of SU{2) - this valence-bond formation is clearly a crucial 
feature determining the structure of the quantum disordered phase. It is well-known 
that this is impossible for SU{2N) for N > 1 - there is no generalization of the 
second-rank, antisymmetric, invariant tensor e to general SU (27V). 

The proper generalization turns out to be to the group Sp{N) p4) . This group is 
defined by the set of 2N x 2N unitary matrices U such that 

U^JU = J (23) 



where 



/ 1 

-1 



(24) 



V 



7 



is the generalization of the e tensor to TV > 1; it has N copies of e along the 
diagonal. It is clear that Sp{N) C SU{2N) for > 1, while Sp{l) ^ SU{2). 
The bf bosons transform as the fundamental representation of Sp{N); the "spins" 
on the lattice therefore belong to the symmetric product of fundamentals, which 
is also an irreducible representation. Valence bonds 



J^^b^- b^. 



(25) 



can be formed between any two sites; this operator is a singlet under Sp{N) because 
of (23 1. The form (22 1 of "H has a natural generalization to general Sp{N): 



i>3 



2N 



J' 



(26) 
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Fig. 9 Phase diagram of the 2D Sp{N) antiferromagnet "H as a function of the "spin" 715 ; from 
Refs. |54||59}j62) The "quantum disordered" region preserves Sp(A'') spin rotation invariance, and 
there is no magnetic long-range order; however, the ground states here have new types of emergent 
order (VBS or Z2 topological order), which are described in the text. On the square lattice, the Z2 
spin liquid phases also break a global lattice rotational symmetry, and so they have 'Ising-nematic' 
order; the Z2 spin liquids on the triangular and kagome lattices do not break any lattice symmetry. 



where the indices a, /?, 7. 6 now run over 1 . . . 2N . We recall also that the constraint 
(17 1 must be imposed on every site of the lattice. 

We now have a two-parameter (rif,, N) family of models T-L for a fixed realization 
of the Jij. It is very instructive to consider the phase diagram of H as a function of 
these two parameters (Fig.[9]l. 

The limit of large n^, with N fixed leads to the semi-classical theory. For the 
special case of SU{2) antiferromagnets with a two-sublattice collinear Neel ground 
state, the semiclassical fluctuations are described by the 0(3) non-linear sigma 
model. For other models ||59 63 -68), the structure of the non-linear sigma mod- 
els is rather more complicated and will not be considered here. 

A second limit in which the problem simplifies is large at fixed rif, |59 69 1. 
It can be shown that in this limit the ground state is quantum disordered. Fur- 
ther, the low-energy dynamics of H is described by an effective quantum-dimer 
model 1 ,59, 70| , with each dimer configuration representing a particular pairing of 
the sites into valence-bonds. There have been extensive studies of such quantum 
dimer models which we will not review here. All the quantum dimer model stud- 
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ies in the "quantum disordered" region of Fig. |9] have yielded phases which were 
obtained earlier p4) by the methods to be described below. 

The most interesting solvable limit is obtained by fixing the ratio of rit and N 

n = - (27) 



and subsequently taking the limit of large N p8) ; this limit will be studied in this 
section in considerable detail. The implementation of H in terms of bosonic opera- 
tors also turns out to be naturally suited for studying this limit. The parameter k is 
arbitrary; tuning k modifies the slope of the line in Fig. [9] along which the large N 
limit is taken. From the previous limits discussed above, one might expect that the 
ground state of T-L has magnetic long range order (LRO) for large k, and is quantum- 
disordered for small n. We will indeed find below that for any set of Jij there is a 
critical value of n — k.c which separates the magnetically ordered and the quantum 
disordered phase. 

The transition at second-order at = oo, and is a powerful feature 

of the present large-A^ limit. In the vicinity of the phase transition, we expect the 
physics to be controlled by long -wavelength, low-energy spin fluctuations; the large- 
TV method offers an unbiased guide in identifying the proper low-energy degress of 
freedom and determines the effective action controlling them. Having obtained a 
long-wavelength continuum theory near the transition, one might hope to analyze 
the continuum theory independently of the large-A^ approximation and obtain re- 
sults that are more generally valid. 

We will discuss the structure of the N — oo mean-field theory , with — kN 



in Section 2.2.1 The long-wavelength effective actions will be derived and used to 



describe general properties of the phases and the phase transitions in Section 2.2.2 



2.2.1 Mean-field theory 

We begin by analyzing T-L si N ^ oo with n\j — kN. As noted above, this limit is 
most conveniently taken using the bosonic operators. We may represent the partition 
function of H by 



VQVbVX exp 



CdT , 



(28) 



where 



E 



N 



(29) 



18 



Subir Sachdev 



Here the fix the boson number of at each site; r-dependence of all fields is im- 
pUcit. The complex field Q was introduced by a Hubbard-Stratonovich decoupling 
of T-L: performing the functional integral over Q reproduces the exchange coupling 



in Eq. (26 1. An important feature of the lagrangian C is its U (1) gauge invariance 



under which 

^ C(«)exp(zp,(T)) 
Qi.j Qij exp (-ipi(r) - ipj{T)) 

\^K + ^{t). (30) 

The functional integral over C faithfully represents the partition function apart from 
an overall factor associated with this gauge redundancy. 

The 1 /N expansion of the free energy can be obtained by integrating out of C 
the 2A^-component h,b fields to leave an effective action for Q, X having coefficient 

(because Ub oc N). Thus the — > oo limit is given by minimizing the effective 
action with respect to "mean-field" values of Q ~ Q, iX — X (we are ignoring here 
the possibility of magnetic LRO which requires an additional condensate = {b") 
- this has been discussed elsewhere |54 62j). This is in turn equivalent to solving 
the mean-field Hamiltonian 

Hmf = E (n'^AI _ ^j^^bfb^ + H.c) 

+ J2~^^{blbt-nb). (31) 

i 

This Hamiltonian is quadratic in the boson operators and all its eigenvalues can be 
determined by a Bogoliubov transformation. This leads in general to an expression 
of the form 

-Hmf = Emf [Q, A] + E [Q, XhU7" ■ (32) 

The index /i extends over 1 . . . number of sites in the system, Em f is the ground 
state energy and is a functional of Q, X, uj^ is the eigenspectrum of excitation ener- 
gies which is a also a function of Q, X, and the 7^ represent the bosonic eigenop- 
erators. The excitation spectrum thus consists of non-interacting spinor bosons. The 
ground state is determined by minimizing Emf with respect to the Qij subject to 
the constraints 

%^=0. (33) 
oXi 

The saddle-point value of the Q satisfies 

Q^J = ( Jc/^^r^f .) (34) 
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Note that Qij — —Qji indicating that Qij is a directed field - an orientation has to 
be chosen on every link. 

We now describe the ground state configurations of the Q, A fields and the nature 
of the bosonic eigenspectrum for the J1-J2-J3 model. We examined the values of 
the energy Emf for Qij configurations which had a translational symmetry with 
two sites per unit cell. For all parameter values configurations with a single site per 
unit cell were always found to be the global minima. We will therefore restrict our 
attention to such configurations. The Aj field is therefore independent of i, while 
there are six independent values of Qij : 





= Qi,x 


Qi,i--\~y 


= Ql,y 


*i,i-\-y-{~x 






— Q2,y — x 


Qi.i-\-2x 


^ Qs.x 


Qi,i^2y 


^ Qs.y ■ 



(35) 

For this choice, the bosonic eigenstates are also eigenstates of momentum with mo- 
menta k extending over the entire first Brillouin zone. The bosonic eigenenergies 
are given by 

Ak = Ji {Qi,x sinfcj; + Qi^y sinfcy) 

+J2 [Q2,y+x su\{ky + fc^) + Q2,y-x sin(fcy - fcj.)) 

+-^^3 (Q3,x sin(2fc^) + Qa^y sin(2fcj^)) . (36) 

We have numerically examined the global minima of Em f as a function of the 
three parameters J2/J1, J3/J1, and N/rif, p4l|62l. The values of the Qij at any 
point in the phase diagram can then be used to classify the distinct classes of states. 



The results are summarized in Figs. 10 and 11 which show two sections of the 



three-dimensional phase diagram. All of the phases are labeled by the wavevector at 
which the spin structure factor has a maximum. This maximum is a delta function 
for the phases with magnetic LRO, while it is simply a smooth function of k for the 
quantum disordered phases (denoted by SRO in Figs. [T0| and [TT]i. The location of 
this maximum will simply be twice the wavevector at which lo^ has a mimimum: 
this is because the structure factor involves the product of two bosonic correlation 
functions, each of which consists of a propagator with energy denominator w^. 

Each of the phases described below has magnetic LRO for large Ub/N and is 
quantum disordered for small ni,/N. The mean-field result for the structure of all 
of the quantum disordered phases is also quite simple: they are featureless spin 
fluids with free spin- 1/2 bosonic excitations ("spinons") with energy dispersion ujk 
which is gapped over the entire Brillouin zone. Some of the quantum disordered 
phases break the lattice rotation symmetry (leading to Tsing-nematic' order) even 
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Fig. 10 Ground states of the J1 — J2 — J3 model for J3 = as afunction of J2/J1 and N/ni, (ni, = 
2S for SU{2)). Thick (thin) lines denote first (second) order transitions at Af = 00. Phases are 
identified by the wavevectors at which they have magnetic long-range-order (LRO) or short-range- 
order (SRO); the SRO phases are "quantum disordered" as in Fig. [9] The links with Qp 7^ in 
each SRO phase are shown. The large N/rii,, large J2/J1 phase has the two sublattices decoupled 
at Af = 00. All LRO phases above have two-sublattice collinear Neel order. All the SRO phases 
above have valence bond solid (VBS) order at finite N for odd rn; this is illustrated by the thick, 
thin and dotted lines. 

at N = 00 - these will be described below. The mininum energy spinons lie at a 
wavevector fco and decreases as rn,/N . The onset of magnetic LRO occurs at 
the value of rn,/N at which the gap first vanishes: Wfe,, = 0. At still larger values of 
rib/N, we get macroscopic bose condensation of the b quanta at the wavevector fco^ 
leading to magnetic LRO at the wavevector 2fco. 

We now turn to a description of the various phases obtained. They can be broadly 
classified into two types; 

Commensurate collinear phases 



In these states the wavevector fco remains pinned at a commensurate point in the 
Brillouin zone, which is independent of the values of J2/J1, J3/J1 and nt,/N. In 
the LRO phase, the spin condensates on the sites are either parallel or anti-parallel 
to each other, which we identify as collinear ordering. This implies that the LRO 
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Fig. 11 As in Fig.[lO| but for J3/ Ji = 0.35. The (0, tt)sro and (tt, tt)sro phases have VBS 
order as illustrated in Fig.^] The {q, q)sRO and (9i t^)sro phases are Z2 spin liquids: they have 
topological order, and a topological 4-fold degeneracy of the ground state on the torus. The Z2 
spin liquids here also have Ising-nematic order, i.e., they break the 90° rotation symmetry of the 
square lattice, which leads to an additional 2-fold degeneracy. The {q, q)LRO and (g, ■k)i^ro have 
magnetic long-range order in the form of an incommensurate spiral. The 2 shaded circles indicate 
regions which map onto the generalized phase diagram in Fig.|12| 



phase remains invariant under rotations about the condensate axis and the rotation 
symmetry is not completely broken. 

Three distinct realizations of such states were found 



a. (tt, tt) 

This is the usual two-sublattice Neel state of the unfrustrated square lattice and its 
quantum-disordered partner. These states have 



}l,x — Ql,y 7^ Oj Q2,y+x — Q2,y-x — Q3,x — Qz,v — • 



(37) 



From (36 1, the minimum spinon excitation occurs at feg = ±(7r/2, 7r/2). The SRO 
states have no broken symmetry at — 00. The boundary between the LRO and 
SRO phases occurs at N/ni, < 2.5, independent of J2/ Ji (Fig. 10 1. This last feature 
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is surely an artifact of the large N limit. Finite N fluctuations should be stronger as 
J2/J1 increases, causing the boundary to bend a little downwards to the right. 



b. (tt, 0) or (0, tt) 

The (0, tt) states have 

Qi,x = 0, Qi^y ^ 0, Q2,y+x = Q2,y-x ^ 0, and Qs^^ = Q^^y = (38) 

and minimum energy spinons at feo = ±(0, 7r/2). The degenerate (7r,0) state is 
obtained with the mapping x ^ y. The SRO state has a two-field degeneracy due 
to the broken x ^ y lattice symmetry: the order associated with this symmetry is 
referred to as 'Ising-nematic' order We can use the Q variables here to define an 
Ising nematic order parameter 

IQI.P- IQlyP. (39) 

This is a gauge-invariant quantity, and the square lattice symmetry of the Hamilto- 
nian implies that (I) = unless the symmetry is spontaneously broken. The sign 
of (I) chooses between the (tt, 0) and (0, tt) states. The LRO state again has two- 
sublattice collinear Neel order, but the assignment of the sublattices is different from 
the (tt, tt) state. The spins are parallel along the x-axis, but anti-parallel along the 
y-axis. 

An interesting feature of the LRO state here is the occurrence of "order-from- 
disorder" |71 1. The classical limit ijib/N = oo) of this model has an accidental de- 
generacy for J2I J\ > 1/2: the ground state has independent collinear Neel order on 
each of the A and B sublattices, with the energy independent of the angle between 
the spins on the two sublattices. Quantum fluctuations are included self-consistently 
in the N — 00, rib/N finite, mean-field theory, and lead to an alignment of the spins 
on the sublattices and LRO at (0, vr). The orientation of the ground state has thus 
been selected by the quantum fluctuations. 

The (0, tt) states are separated from the (tt, tt) states by a first-order transition. In 
particular, the spin stiffnesses of both states remain finite at the boundary between 
them. This should be distinguished from the classical limit in which the stiffness 
of both states vanish at their boundary J2 — Ji/2; the finite spin stiffnesses are 
thus another manifestation of order-from-disorder. At a point well away from the 



singular point J2 = Ji/2, rih/N = cx) in Fig. 10 the stiffness of both states is 
of order N{ni,/N)'^ for N = 00 and large rih/N; near this singular point however 
the stiffness is of order N{ni,/N) and is induced purely by quantum fluctuations. 
These results have also been obtained by a careful resummation of the semiclassical 
expansion |[72l|73|. 
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c. "Decoupled" 



For Jij J\ and iV/nt both large, we have a "decoupled" state (Fig. 10 1 with 

Q2,y^x = Q2,y~x + and Qx^Qi = 0. (40) 

In this case Qp is non-zero only between sites on the same sublattice. The two 
sublattices have Neel type SRO which will be coupled by finite N fluctuations. The 
N — oo state does not break any lattice symmetry. This state has no LRO partner. 



Incommensurate phases 

In these phases the wavevector feo and the location of the maximum in the structure 
factor move continuously with the parameters. The spin-condensate rotates with a 
period which is not commensurate with the underlying lattice spacing. Further the 
spin condensate is coplanar. the spins rotate within a given plane in spin space and 
are not coUinear. There is no spin rotation axis about which the LRO state remains 
invariant. 

Further, no states in which the spin condensate was fully three dimensional 
("double-spiral" or chiral states) were found; these would be associated with com- 
plex values of Qp. All the saddle points possesed a gauge in which all the Qp were 
real. Time-reversal symmetry was therefore always preserved in all the SRO phases 
ofFigs.[TO]and[TT] 

The incommensurate phases occur only in models with a finite J3 (Fig. Ill, at 
least at = 00. There were two realizations: 



d. (tt,?) or (q,7r) 

Here q denotes a wavevector which varies continuously between and tt as the 
parameters are changed. The (g, tt) state has 

Qi,x ^ Qi,y 7^ 0, Q2,x+y = Q2,y-x ^ 0, Qz.x ^ and Qa^j^ = 0; (41) 

the degenerate (tt, q) helix is obtained by the mapping x ^ y. The SRO state has 
a two-fold degeneracy due to the broken x ^ y lattice symmetry, and so this state 



has Ising-nematic order. The order parameter in Eq. ( 39 1 continues to measure this 
broken symmetry. 



e. {q, q) or (q, -q) 

The ((7, q) state has 

Qi,x = Qi,y ^ 0, Q2,x+y 0, Q2^y-x = 0, Qs,^ = Qa.y ^ 0; (42) 
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this is degenerate with the {q, —q) phase. The SRO state therefore has a two-fold 
degeneracy due to a broken lattice reflection symmetry, and so it also has Ising 
nematic order. However, the Ising symmetry now corresponds to reflections about 
the principle square axes, and the analog of Eq. ( [39| is now 

^^\Q2,x+y\^ -\Q2,y-x\^- (43) 

As we noted above, the broken discrete symmetries in states with SRO at (0, tt) 
and {q, tt) are identical: both are two-fold degenerate due to a breaking of the x ^ y 
symmetry. The states are only distinguished by a non-zero value of in the {q, tt) 
phase and the accompanying incommensurate correlations in the spin-spin correla- 
tion functions. However Q3 is gauge-dependent and so somewhat unphysical as an 
order parameter. In the absence of any further fluctuation-driven lattice symmetry 
breaking, the transition between SRO at (0,7r) and (g, tt) is an example of a dis- 
order line |74J; these are lines at which incommensurate correlations first turn on. 
However, we will see that quantum fluctuations clearly distinguish these two phases, 
which have confined and deconfined spinons respectively, and the associated topo- 
logical order requires a phase transitions between them. 

An interesting feature of Fig. [TT] is that the commensurate states squeeze out 
the incommensurate phases as N /uf, increases. We expect that this suppression of 
incommensurate order by quantum fluctuations is a general feature of frustrated 
antiferromagnets. 



2.2.2 Fluctuations - long wavelength effective actions 



We now extend the analysis of Section 2.2.1 beyond the mean-field theory and ex 



amine the consequences of corrections at finite N. The main question we hope to 
address are: 

• The mean-field theory yielded an excitation spectrum consisting of free spin- 1/2 
bosonic spinons. We now want to understand the nature of the forces between 
these spinons and whether they can lead to confinement of half-integer spin ex- 
citations. 

• Are there any collective excitations and does their dynamics modify in any way 
the nature of the mean field ground state ? 

The structure of the fluctuations will clearly be determined by the low-energy 
excitations about the mean-field state. We have already identified one set of such 
excitations: spinons at momenta near mimima in their dispersion spectrum, close to 
the onset of the magnetic LRO phase whence the spinon gap vanishes. An additional 
set of low-lying spinless excitations can arise from the fluctuations of the Qij and 



fields about their mean-field values. The gauge-invariance ( 30 1 will act as a powerful 
restriction on the allowed terms in the effective action for these spinless fields. We 
anticipate that the only such low-lying excitations are associated with the and the 
phases of the Qij. We therefore parametrize 
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QiA+e^ = QiA+ip exp i~iOp) , (44) 

where the vector ip connects the two sites of the lattice under consideration, Q is 
the mean-field value, and Op is a real phase. The gauge invariance ( 30 1 implies that 
the effective action for the Op must be invariant under 

Op^Op+p^+ . (45) 

Upon performing a Fourier transform, with the link variables Op placed on the cen- 
ter of the links, the gauge invariance takes the form 

Op(k) -5> Opik) + 2p(k) cos(fcp/2) , (46) 

where kp = k ■ Bp. This invariance implies that the effective action for the Op, after 
integrating out the b quanta, can only be a function of the following gauge-invariant 
combinations: 



Ipq = 2cos{kq/2)Op{k) -2cos{kp/2)Oq{k) . (47) 

We now wish to take the continuum limit at points in the Brillouin zone where the 
action involves only gradients of the Op fields and thus has the possibility of gapless 
excitations. This involves expanding about points in the Brillouin zone where 

cos(/cp/2) = for the largest numbers of Cp. (48) 

We now apply this general principle to the J1-J2-J3 model. 



Commensurate colUnear phases 



We begin by examining the (tt, 7r)-SRO phase. As noted in ( 37 1, this phase has the 



mean field values Qi^x — Qi.y 7^ 0, and all other Qij zero. Thus we need only 



examine the condition (48 1 with ip = ex,ey. This uniquely identifies the point 



k = G = (tt, tt) in the Brillouin zone. We therefore parametrize 

©,(r) = e^«M,(r) (49) 

and similarly for Oy\ it can be verified that both O and are real in the above 
equation. We will also be examining invariances of the theory under gauge transfor- 
mations near G: so we write 

p{r) = e^« '-C(r) . (50) 



It is now straightforward to verify that the gauge transformations ( |46] l are equivalent 
to 

Ax^Ax + dxC (51) 
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and similarly for Ay. We will also need in the continuum limit the component of A 
near the wavevector G. We therefore write 



iX, = X + ie"^-''Ar{ri). 
Under gauge transformations we have 



(52) 



(53) 



Thus Ax, Ay, Ar transform as components of a continuum U{V) vector gauge field. 

We will also need the properties of the boson operators under the gauge transfor- 
mation C^. From pO] ) and ( [50| we see that the bosons on the two sublattices (A, E) 
carry opposite charges ±1: 



hA -> hAe 



Finally, we note that the bosonic eigenspectrum has a minimum near k 
(7r/2, tt/2); we therefore parametrize 



(54) 
fco = 

(55) 



We insert the continuum parameterizations ( |49| l, ( 52 1 and ( 55 i into the functional 
integral ([29]l, perform a gradient expansion, and transform the Lagrangian C into 



C = 



dT 



- iAr V'2a 



+ A d^n' + 1^-201') - 4JiQi (V'?^2a + rici^T) 

+JiQia^ [(V + iA) (V - lA) V'2a 



(56) 



We now introduce the fields 



= (^?+^r)/V2 

^" = «-V2"*)/V2. 

Following the definitions of the underlying spin operators, it is not difficult to show 
that the Neel order parameter ipa (which is proportional to n in Q) is related to the 

Za by 

= zy^'^pz^ . (57) 



From Eq. (56 1, it is clear that the the tt fields turn out to have mass A + AJiQi, 



while the z fields have a mass A— 4Ji(5i which vanishes at the transition to the LRO 
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phase. The tt fields can therefore be safely integrated out, and H yields the following 
effective action, valid at distances much larger than the lattice spacing 1 60 61]: 



^cff^y^y^ dr|i(9^-*A^)z"|2 + — |z"in . (58) 

Here ^ extends over x, y,T, c — ^/SJiQia is the spin-wave velocity, A — (A^ — 
IGJfQiY^^ is the gap towards spinon excitations, and Af — Arjc. Thus, in its 
final form, the long-wavelength theory consists of a massive, spin-1/2, relativistic, 
boson (spinon) coupled to a compact C/(l) gauge field. By 'compact' we mean 
that values and + 27r are identified with each other, and the gauge field lives 



on a circle: this is clearly required by Eq. (44 1 



At distances larger than cj A, we may safely integrate out the massive z quanta 
and obtain a a compact gauge theory in 2h-1 dimensions. This theory was 
argued by Polyakov | |75|[76) to be permanently in a confining phase, with the con- 
finement driven by "monopole" tunnelling events. The compact J7(l) gauge force 
will therefore confine the z" quanta in pairs. So the conclusion is that the (tt, tt) g^o 
does not possess 5=1/2 spinon excitations, as was the case in the mean field the- 
ory. Instead, the lowest-lying excitations with non-zero spin will be triplons, similar 
to those in Section [2T| A further important effect here, not present in the [/(I) 
gauge theories considered by Polyakov, is that the monopole tunnelling events carry 
Berry phases. The influence of these Berry phases has been described ||60j[6l | and 
reviewed |39 | elsewhere, and so will not be explained here. The result is that the con- 
densation of monopoles with Berry phases leads to valence bond solid (VBS) order 
in the ground state. This order is associated with the breaking of the square lattice 
space group symmetry, as illustrated in Figs. 10 and 12 below. For the {'k^'k)sro 



phase, this means that the singlet spin correlations have a structure similar to that in 
Fig. |4] In other words, the square lattice antiferromagnet spontaneously acquires a 
ground state with a symmetry similar to that of the paramagnetic phase of coupled- 
dimer antiferromagnet. Because the VBS order is spontaneous, the ground state is 
4-fold degenerate (associated with 90° rotations about a lattice site), unlike the non- 
degenerate ground state of the dimerized antiferromagnet of Section [2T| VBS states 
with a plaquette ordering pattern can also appear, but are not shown in the figures. 
The quantum phase transition between the (tt, 'k)sro and (tt, 'n^LRO phases has 



been the topic of extensive study. The proposal of Refs. |77 78] is that monopoles 
are suppressed precisely at the quantum critical point, and so the continuum action 
in Eq. ( [58] l constitutes a complete description of the critical degrees of freedom. 
It has to be supplemented by a quartic non-linearity (|z"P) , because such short- 
range interactions are relevant perturbations at the critical point. A review of this 
deconfined criticality proposal is found elsewhere p7) . 

The properties of the (0, vr) phase are very similar to those of the (tt, vr) phase 
considered above. It can be shown quite generally that any quantum disordered state 
which has appreciable commensurate, collinear spin correlations will have similar 
properties: confined spinons, a collective mode described by a compact UiV) gauge 
field, and VBS order for odd rib. 
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Incommensurate phases 



We now turn to a study of the incommensurate phases. It is not difficuh to show 



that in this case it is not possible to satisfy the constraints (48 1 at any point in the 
Brillouin zone for all the non-zero Qp. This implies that, unlike the commensurate 
phases, there is no gapless collective gauge mode in the gaussian fluctuations of 
the incommensurate SRO phases. This has the important implication that the mean- 
field theory is stable: the structure of the mean-field ground state, and its spinon 
excitations will survive fluctuation corrections. Thus we obtain a stable 'spin liquid' 
with bosonic S = 1/2 spinon excitations. We will now show that these spinons carry 
a Z2 gauge charge, and so this phase is referred to as a Z2 spin liquid. The Z2 gauge 
field also accounts for 'topological order' and a 4-fold ground state degeneracy on 
the torus. 

The structure of the theory is simplest in the vicinity of a transition to a com- 
mensurate collinear phase: we now examine the effective action as one moves from 



the (tt, 7r)-SRO phase into the {q, g)-SRO phase (Fig. 1 1 , a very similar analysis can 
be performed at the boundary between the (vr, 7r)-SRO and the (tt, (7)-SR0 phases). 
This transition is characterized by a continuous turning on of non-zero values of 



Qi.i+2x and Qi_i+2y- It is easy to see from Eq. ( 30 1 that these fields trans- 
form as scalars of charge ±2 under the gauge transformation associated with A^. 
Performing a gradient expansion upon the bosonic fields coupled to these scalars we 
find that the Lagrangian C of the (tt, 7r)-SR0 phase gets modified to 

(Aa ■ (ja/sV-rVVf) + Ab ■ (j"''V2aV^2/3) +C.C.) , (59) 

where Aa,b ai'e two-component scalars ee {JzQ'a^x + J2Q2,y+x, JaQs.y + J2 
Q2,y+x) with the sites on the ends of the link variables on sublattices A^B. Fi- 
nally, as before, we transform to the z, tt variables, integrate out the tt fluctuations 
and obtain l|62l 



'S'off — 



/ Wa C "^^{'^^^ " + + ^ • {Jc^z'^^z^) + c.c. 

+ KA\{d^ + 2iA^,)A\^ + saA^ + terms quartic in z", A I . (60) 



Here Sz = A'^/c?, A = {Aa + A*^) / {2JiQia) is a complex scalar of charge 
—2, and Ka is a stiffness. We have explicitly written the quadratic terms in the 
effective action for the A: these are generated by short wavelength fluctuations of 
the 6" quanta. We have omitted quartic and higher order terms which are needed to 
stabilize the theory when the 'masses' or sa are negative, and are also important 
near the quantum phase transitions. This effective action is also the simplest theory 
that can be written down which couples a spin-1/2, charge 1, boson z", a compact 
U{1) gauge field A^, and a two spatial component, charge —2, spinless boson A. 
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7.2 spin liquid 

(^")=0 , (A>^0 
Coplanar SRO 



, (A)/0 
Coplanar LRO 




V6S 



(^«) = , (A) = 
Collinear SRO 







(^")/0 , (A) = 
Collinear LRO 
\V \\ 



5a 



Fig. 12 Phase diagram of the theory 5off in Eq. |60| for the bosonic spinons and the charge 
-2 spinless boson A. Fig.[TO]contains examples of the region Sj\ > 0. Fig.[TT]contains 2 separate 
instances of 4 phases meeting at a point as above, with the 4 phases falling into the classes labeled 
above; these points are labeled in both figures by the shaded circles. 



It is the main result of this section and summarizes essentially all of the physics we 
are trying to describe. 



We now describe the various phases of Scs, which are summarized in Fig. 12 



1. Commensurate, coHineai", LRO: (z") 0, (A) = 

This is the (tt, tt)lr.o state with commensurate, collinear, magnetic LRO. 

2. Commensurate, collineai', SRO: (z") = 0, (A) = 

This is the (tt, tt)sro "quantum-disordered" state with collinear spin correlations 
peaked at (tTjTt). Its properties were described at length above. The compact 
U{1) gauge force confines the z" quanta. The spinless collective mode associated 
with the gauge fluctuations acquires a gap from monopole condensation, and the 
monopole Berry phases induce VBS order for odd n^. 

3. Incommensurate, coplanar, SRO: (z") — 0, (A) ^ 

This is the incommensurate phase with SRO at (q, q) which we want to study. 
It is easy to see that condensation of A necessarily implies the appearance of 
incommensurate SRO: ignore fluctuations of A about (A) and diagonalize the 
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quadratic form controlling the fluctuations; the minimum of the dispersion of 
the quanta is at a non-zero wavevector 

ko^[{A^),{Ay))/2. (61) 

The spin structure factor will therefore have a maximum at an incommensurate 
wavevector This phase also has a broken lattice rotation symmetry due to the 
choice of orientation in the x — y plane made by A condensate, i.e., it has Ising- 
nematic order. 

The condensation of A also has a dramatic impact on the nature of the force be- 
tween the massive z" quanta. Detailed arguments have been presented by Frad- 
kin and Shenker | ,79J that the condensation of a doubly charged Higgs scalar 
quenches the confining compact U{1) gauge force in 2+\ dimensions between 
singly charged particles. We can see this here from Eq. ( [60| by noticing that the 
condensation of A expels A^ by the Meissner effect: consequently, monopoles in 
A^ are connected by a flux tube whose action grows linearly with the separation 
between monopoles. The monopoles are therefore confined, and are unable to in- 



duce the confinement of the z" quanta. From Eq. ( 60 1 we also see that once A is 
condensed, the resulting theory for the spinous only has an effective Z2 gauge in- 
variance: Ref. |79 1 argued that there is an effective description of this free spinon 
phase in terms of a Z2 gauge theory. The excitation structure is therefore very 
similar to that of the mean-field theory: spin- 1/2, massive bosonic spinous and 
spinless collective modes which have a gap. The collective mode gap is present 
in this case even at iV = 00 and is associated with the condensation of A. 



This state is also ' topologically ordered'. We can see this by noticing |54 55] 



that it carries stable point-like excitations which are 27r vortices in either of A^ 
or Ay. Because of the screening by the gauge field, the vortices carry a finite 
energy (this is analogous to the screening of supercurrents by the magnetic field 
around an Abrikosov vortex in a superconductor). Because the A carry charge 
—2, the total flux trapped by a vortex is tt. Thus the vortices are also stable to 
monopole tunneling events, which change the flux by integer multiples of 2tt. 
A Za spinon circumnavigating such a vortex would pick up an Aharanov-Bohm 
phase factor of tt (because the spinous have unit charge), and this is equivalent 
to the statement that the vortex and the spinon obey mutual 'semionic' statistics. 
All these characteristics identify the vortex excitation as one dubbed later [80] as 
the vison. 

The vison also allows us to see the degeneracy of the gapped ground state on sur- 
faces of non-trivial topology. We can insert a vison through any of the 'holes' in 
surface, and obtain a new candidate eigenstate. This eigenstate has an energy es- 
sentially degenerate with the ground state because the core of the vortex is within 
the hole, and so costs no energy. The 'far field' of the vison is within the system, 
but it costs negligible energy because the currents have been fully screened by 
A^ in this region. Thus we obtain a factor of 2 increase in the degeneracy for 
every 'hole' in the surface (which is in turn related to the genus of the surface). 
The state so obtained is now referred to as a Z2 spin liquid, and has been la- 
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beled as such in the figures. As we noted above, the present theory only yields 
Z2 spin liquids with Ising-nematic order, associated with broken symmetry of 
90° lattice rotations. We also note an elegant exactly solvable model described 
by Kitaev |81 1, which has spinon and vison excitations with the characteristics 
described above, but without the Ising-nematic order. 
4. Incommensurate, coplanar, LRO: (z") 7^ 0, (A) ^ 

The condensation of the z quanta at the wavevector feo above leads to incom- 
mensurate LRO in the (g, q)hRO phase, with the spin condensate spiraling in the 
plane. 

We also note a recent work p7][82] [83) which has given a dual perspective on the 
above phases, including an efficient description of the phase transitions between 
them, and applied the results to experiments on k-(ET)2Cu2(CN)3. 



3 d-wave superconductors 

In our discussion of phase transitions in insulators we found that the low-energy 
excitations near the critical point were linked in some way to the broken symmetry 
of the magnetically ordered state. In the models of Section[2T|the low-energy exci- 



tations involved long wavelength fluctuations of the order parameter In Section 2.2 



the connection to the order parameter was more subtle but nevertheless present: the 



field Za in Eq. ( 60 1 is a 'fraction' of the order parameter as indicated in ( 57 1, and the 
gauge field represents a non-coplanarity in the local order parameter orientation. 

We will now move from insulators to the corresponding transitions in d-wave su- 
perconductors. Thus we will directly address the criticality of the magnetic QPT at 
X = Xs VS\ Fig.[T] We will also consider the criticality of the 'remnant' Ising-nematic 
ordering ■aS.x — Xm within the superconducting phase. A crucial property of d-wave 
superconductors is that they generically contain gapless, fermionic Bogoliubov ex- 
citations, as we will review below. These gapless excitations have a massless Dirac 
spectrum near isolated points in Brillouin zone. While these fermionic excitations 
are present in the non-critical d-wave superconductor, it is natural to ask whether 
they modify the theory of the QPT. Even though they may not be directly related 
to the order parameter, we can ask if the order parameter and fermionic excitations 
couple in interesting ways, and whether this coupling modifies the universality class 
of the transition. These questions will be answered in the following subsections. 

We note that symmetry breaking transitions in graphene are also described by 
field theories similar to those discussed in this section I 



3.1 Dirac fermions 



We begin with a review of the standard BCS mean-field theory for a d-wave su- 
perconductor on the square lattice, with an eye towards identifying the fermionic 
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Bogoliubov quasiparticle excitations. For now, we assume we are far from any QPT 
associated with SDW, Ising-nematic, or other broken symmetries. We consider the 
generaHzed Hamiltonian 



tj 



= ^efc4aCfc" + Ji'^Si ■ Sj 



(62) 



where Cja is the annihilation operator for an electron on site j with spin a =t, i, 
Cka is its Fourier transform to momentum space, Sk is the dispersion of the electrons 
(it is conventional to choose Sk — ~2ti{cos{kx) + cos{ky)) — 2t2(cos(fca; + ky) + 
cos{kx — ky)) — jjL, with <i 2 the first/second neighbor hopping and /i the chemical 



potential), and the Ji term is the same as that in Eq. ( 15 1 with 



Sja 



-It 



(63) 



and cr" the Pauli matrices. We will consider the consequences of the further neighbor 



exchange interactions in ( 15 1 for the superconductor in Section 3.3.1 below. Apply- 
ing the BCS mean-field decoupling to Htj we obtain the Bogoliubov Hamiltonian 



Hbcs 



k 



h.c. 



(64) 



For a wide range of parameters, the ground state energy is optimized by a d^2_y2 
wavefunction for the Cooper pairs: this corresponds to the choice ~ — 
A^2_y2. The value of Ax2_y2 is determined by minimizing the energy of the BCS 
state 

d^k 



EbCS = Jl\^x^-t, 



47r2 



[Ek - ek] 



(65) 



where the fermionic quasiparticle dispersion is 



Ek 



el 



\JiAx2_^y2(coskx — cos ky)\ 



1/2 



(66) 



The energy of the quasiparticles, Ek, vanishes at the four points (±Q, ±Q) at 
which Ek = 0. We are especially interested in the low-energy quasiparticles in the 
vicinity of these points, and so we perform a gradient expansion of Hbcs near each 
of them. We label the points Qi = (Q,Q), Q2 = {-Q,Q), Q3 = {-Q,-Q), 
Qi = {Q, ^Q) and write 



iQ2rj 



+ /4a(r,)e'«-'■^ 



(67) 



while assuming that the /1-4. c((r) are slowly varying functions of r. We also intro- 
duce the bispinors !fi = (/i^, /]^, /i^, ~fl^), and 'P2 = {f2t, fl^, hi, -fl-^)^ and 
then express Hbcs in terms of S'i.2 while performing a spatial gradient expansion. 
This yields the following effective action for the fermionic quasiparticles: 
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— drd r 



. VF 



Wl dr - i^id, + dy)T' - i^i-d, + dy)r- 



V2' ^ ' " V2 



.(68) 



where the r^'^ are 4 x 4 matrices which are block diagonal, the blocks consisting 
of 2 X 2 Pauli matrices. The velocities vf,a are given by the conical structure of 
near the Qi-^. we have vp = \^ k£k\k=Qa \ and va = |«/i^3;2_j,2\/2sin((5)|. In 
this limit, the energy of the S'l fermionic excitations is — (vp(kx + ky)'^/2 + 



ky)^/2)^^^ (and similarly for ^2), which is the spectrum of massless Dirac 



fermions. 



3.2 Magnetic ordering 

We now focus attention on the QPT involving loss of magnetic ordering within the 
d-wave superconductor at a; = Xs in Fig.[T| As in Section |2] we have to now con- 
sider the fluctuations of the SDW order parameter. We discussed two routes to such 
a magnetic ordering transition in Section |2j one involving the vector SDW order 



parameter in Section 2.1 and the other involving the spinor Za in Section 2.2 In 
principle, both routes also have to be considered in the d-wave superconductor. The 
choice between the two routes involves subtle questions on the nature of fractional- 
ized excitations at intermediate scales which we will not explore further here. These 
questions were thoroughly addressed in Ref. |28| in the context of simple toy mod- 
els: it was found that either route could apply, and the choice depended sensitively 
on microscopic details. In particular, it was found that among the fates of the non- 
magnetic superconductor was that it acquired VBS or Ising-nematic ordering, as 



was found in the models explored in Section 2.2 This is part of the motivation for 
the expectation of such ordering in the regime < x < x„i, as indicated in Fig.[T] 
In the interests of brevity and simplicity, we will limit our discussion of the SDW 



ordering transition here to the vector formulation analogous to that in Section 2. 1 



We have full square lattice symmetry, and so allow for incommensurate SDW or- 



dering similar to the {q, q),{<l, —<l) and (tt, q),{q, tt) states of Section 2.2 Because 



there are two distinct but degenerate ordering wavevectors, the complex order pa- 



rameter in Eq. (14 1 is now replaced by two complex order parameters (l>xa and 
<Pya for orderings along (tt, q) and {q, tt) (the orderings along {q, ±q) can be treated 
similarly and we will not describe it explicitly). These order parameters are related 
to the spin operator by 

Sa{r) = •P.^e'-^-"' + ^yac^'^y'- + c.c. (69) 



where — (q,'^) and Ky = {TT,q). As discussed below Eq. (14 1, depending 
upon the structure of the complex numbers <Pxa, '^ya, the SDW ordering can be 
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Tx 


Ty 


R 


/ 


r 




e"l<Pxa 




■Pya 






'Pya 
^la 


-"Pya 




xa 


ya 

'P20C 




1'2c 










-ryi'2cc 



Table 1 Transformations of the fields under operations wliicli generate the symmetry group: 
Tx,y = translation by a lattice spacing in the x, y directions, R = rotation about a lattice site 
by 90°, / = reflection about the y axis on a lattice site, and T = time reversal. The theory is also 
invariant under spin rotations, with i a vector index and a, 13 spinor indices. We define T as an 
invariance of the imaginary tiine path integral, in which $ ^ 2i transform as the complex conjugates 

of ^i,2i, while 2^ are viewed as independent coinplex Grassman fields which transform as 

1,2a 1,2a 



either collinear (i.e., stripe-like) or spiral. Also, as in Eqs. (39 1 and (43 1, we can 
use these SDW order parameters to also define a subsidiary Ising-nematic order 
parameter 



1= <P. 



(70) 



to measure the breaking of x ^ y symmetry. 

Symmetry considerations will play an important role in our analysis of the 'P^.ya 
order parameters and their coupling to the Dirac fermions. In Table [T] we therefore 
present a table of transformations under important operations of the square lattice 



space group: these are easily deduced from the representations in Eq. ([67|) and ( 69 1 



The effective action for the SDW order parameters has a direct generalization 
from ( [l4| : it can be obtained by requiring invariance under the transformations in 
Table irrand has many more allowed quartic nonlinearities EO): 



(Prdr 



+ \dr'Pya\ + 4 \dy^ya\ + 4 \d^^ya\ + S l^.al' + \^cca 



Ul 

2 



+ {\<l>ya\'y 



U2 

2 



{1>.af + i'Pya) 



■Wi\^a:a\ \^ya\ +W2\^xa^ya\ +W3|^xa<?^a| 



(71) 



Remarkably, a fairly complete 5-loop renormalization group analysis of this model 



has been carried out by De Prato et al. |86|, and reliable information on its critical 
properties is now available. 



(We note parenthetically that Eq. (71 1 concerns the theory of the transition at Xg 
from an SDW ordered state to a d-wave superconductor with the full symmetry of 
the square lattice. However, as we have discussed in Section[T]and in the beginning 
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Fig. 13 The filled circles indicate the positions of the gapless Dirac fennions in the square lattice 
Brillouin zone: these are at wavevectors Qi,2.3,4- An SDW fluctuation scatters a femiion at one of 
the nodes by wavevector to a generic point in the Brillouin zone. The final state of the fermion 
has a high energy, and so such processes are suppressed. 



of Section [3] there could be Ising nematic order in the regime Xg < x < Xm - In this 
case one of <Pxa or (l>ya orderings would be preferred, and we need only consider 
the critical fluctuations of this preferred component. The resulting action for this 
preferred component would then be identical to Eq. ( 14 1, with critical properties as 
inRef. I^J.) 

Now we turn to the crucial issue of the coupling between the 'l'x,ya order pa- 
rameter degrees of freedom in and the massless Dirac fermions in Eq. (68 1. 
Again a great deal follows purely from symmetry considerations. The simplest pos- 
sible terms are cubic 'Yukawa' interaction terms like 'l'xa}Pi^2 etc. However, these 
are generically forbidden by translational invariance, or equivalently, momentum 
conservation. In particular, the transformation of the 'Px.ya under translation by one 



lattice spacing follows from ( [69] l, while those of the if'i 2 follow from (67 1. Unless 
the SDW ordering wavevectors K^.y and the positions of the Dirac nodes Qi,2.3.4 
satisfy certain commensurability conditions, the Yukawa coupling will not be in- 



variant under this translation operation. This is illustrated schematically in Fig. 13 



The observed values of the wavevectors are not commensurate, and so we can safely 
neglect the Yukawa term. 

The absence of the Yukawa coupling suggests that the fixed point theory describ- 
ing the QPT at X = Xg in the superconductor may be S<p in Eq. (71 1 alone, i.e., the 
transition is in the same universality class as the insulator However, to ensure this, 
we have to examine the influence of higher terms coupling the degrees of freedom 
of S(f, and Sq,. The simplest couplings not prohibited by translational invariance are 
associated with operators which are close to net zero momentum in both sectors. 
These are further constrained by the other square lattice space group operations 
in Table [T[ requiring invariance under them shows that the simplest allowed terms 
arellSTl 
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S2=^2 j dTffr (|<^,,p - \1>ya?) {^y^Pi + ■ (72) 

The first term is a fairly obvious 'density-density' coupling between the energies of 
the two systems. The second is more interesting: it involves the Ising nematic order 



X, as measured in the order parameter sector by ( 70 1, and in the fermion sector by 
the bilinear shown above. 

Now we can ask if the fixed point described by the decoupled theory S4, + 
is stable under the perturbations in Si and ^2. This involves a computation of the 
scaling dimensions of the couplings 2 at the decoupled theory fixed point. These 
scaling dimensions were computed to 5 -loop order in Ref. | |87| , and it was found that 
dim[i?i] « —1.0, and dim['i92] ~ —0.1. Thus both couplings are irrelevant, and we 
can indeed finally conclude that the SDW onset transition is described by the same 
theory as in the insulator However, the scaling dimension of the i?2 coupling is quite 
small, indicating that it will lead to appreciable effects. Thus we have demonstrated 
quite generally that it is the Ising-nematic order I which is most efficient in coupling 
the SDW order parameter fluctuations to the Dirac fermions. Note that I was not 
chosen by hand, but was selected by the theory among all other possible composite 
orders of the SDW field 'Px,ya- The near zero scaling dimension of i?2 implies that 
it will induce a Unewidth ~ T in the Dirac fermion spectrum. Moreover, because 
this broadening is mediated by I, the broadening will be strongly anisotropic in 
space | ,88J . 



3.3 Ising transitions 

Now we turn our attention to the vicinity of the point x,n in Fig. [T| Although Xm 
was defined in terms of the SDW transition in the metal at high magnetic fields, we 
have also argued in Section [T] and in the beginning of Section [3] that there can also 
be transitions associated with VBS or Ising nematic order near Xm but within the 
superconducting phase at zero field. Strong evidence for a nematic order transition 
near Xm has emerged in recent experiments | [T2][30l[3T[ . 

This section will therefore consider the theory of Ising-nematic ordering within a 
d-wave superconductor Unlike the situation in Section [l!2] we will find here that the 
order parameter and the Dirac fermions are strongly coupled, and the universality 
class of the transition is completely changed by the presence of the Dirac fermions. 



In Section 3.2 we found that although the fermions were moderately strongly cou- 
pled to the critical theory, they were ultimately reduced to spectators to the asymp- 
totic critical behavior. 

Before considering the Ising nematic transition, we will take a short detour in 



Section 3.3.1 and describe another Ising transition associated with the breaking of 



time-reversal symmetry in a d-wave superconductor. This leads to a model which 
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has a somewhat simpler structure, and for which conventional renormahzation group 



techniques work easily. We will return to Ising-nematic ordering in Section 3.3.2 



3.3.1 Time-reversal symmetry breaking 

We will consider a simple model in which the pairing symmetry of the supercon- 
ductor changes from dx2_y2 to d^2_y2 ± id^y. The choice of the phase between the 
two pairing components leads to a breaking of time-reversal symmetry. Studies of 
this transition were originally motivated by the cuprate phenomenology, but we will 
not explore this experimental connection here because the evidence has remained 
sparse. 

The mean field theory of this transition can be explored entirely within the con- 
text of BCS theory, as we will review below. However, fluctuations about the BCS 
theory are strong, and lead to non-trivial critical behavior involving both the col- 
lective order parameter and the Bogoliubov fermions: this is probably the earliest 
known example [89 90| of the failure of BCS theory in two (or higher) dimensions 
in a superconducting ground state. At T > 0, this failure broadens into the "quantum 
critical" region. 



We extend Htj in Eq. (62 1 so that BCS mean-field theory permits a region with 
dxy superconductivity. It turns out that the frustrating interactions as in Eq. ( 15 i are 
precisely those needed. With a J2 interaction, Eq. (|62]) is modified to: 



k (ij) nnn ij 

We will follow the evolution of the ground state of Htj as a function of J2/ Ji- 



The mean-field Hamiltonian is now modified from Eq. (64 1 to 



2 



2 

3,1^ 



where the second summation over v is along the diagonal neighbors x + y and 
—x + y. To obtain d^y pairing along the diagonals, we choose A^+y — —A-^+y = 
Axy We summarize our choices for the spatial structure of the pairing amplitudes 



(which determine the Cooper pair wavefunction) in Fig. 14 The values of Ar^2_y2 
and A^y are to be determined by minimizing the ground state energy (generalizing 



Eq. (65 I) 



d^k 

EbCS ^ .h\Ax^^y^\'' + J2\Axy\'^ - / ^ [^fc - ^fe] , (75) 



where the quasiparticle dispersion is now (generalizing Eq. (66 1) 
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£j, + \JiAx2_y2 (cos kx — cos ky) + 2J2Axy sin kx sin ky 



nl/2 



(76) 



Notice that the energy depends upon the relative phase of Ax2^y2 and A^y'- this 
phase is therefore an observable property of the ground state. 

It is a simple matter to numerically carry out the minimization of Eq. (|76]l, and the 
results for a typical choice of parameters are shown in Fig. 15 as a function J2/J1. 



One of the two amplitudes Ax2^y2 or A^y is always non-zero and so the ground 
state is always superconducting. The transition from pure (1^2 ^y2 superconductivity 
to pure dxy superconductivity occurs via an intermediate phase in which both order 
parameters are non-zero. Furthermore, in this regime, their relative phase is found 
to be pinned to ±7r/2, i.e. 

aig(Axy) = arg(Z\j.2_y2) ± 7r/2 . (77) 



The reason for this pinning can be intuitively seen from Eq. (76i: only for these 
values of the relative phase does the equation — Q never have a solution. In 
other words, the gapless nodal quasiparticles of the (1^2 -y2 superconductor acquire 
a finite energy gap when a secondary pairing with relative phase ±7r/2 develops. By 
a level repulsion picture, we can expect that gapping out the low-energy excitations 
should help lower the energy of the ground state. The intermediate phase obeying 
Eq. (|77jl is called a dx2^y2 + id^y superconductor. 



The choice of the sign in Eq. ( 77 1 leads to an overall two-fold degeneracy in the 
choice of the wavefunction for the dx2-y2 + idxy superconductor. This choice is re- 
lated to the breaking of time-reversal symmetry, and implies that the dx2-y2 + idxy 
phase is characterized by the non-zero expectation value of a Z2 Ising order param- 
eter; the expectation value of this order vanishes in the two phases (the (i^2_y2 and 
dxy superconductors) on either side of the dx2-y2 + idxy superconductor As is con- 
ventional, we will represent the Ising order by a real scalar field 0. Fluctuations of (f> 
become critical near both of the phase boundaries in Fig. [15] As we will explain be- 
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Fig. 15 BCS solution of the phenomenological Hamiltonian Htj in Eq. 173). Shown are the op- 



timum values of the pairing amplitudes |zi^2_„2| and \Axy\ as a function of J2 for ti 



1, 



t2 = —0.25, fi = —1.25, and Ji fixed at Ji = 0.4. The relative phase of the pairing amplitudes 
was always found to obey Eq. ^77) . The dashed lines denote locations of phase transitions between 
d^2_y2, d^2_y2 +idxy, and d^y superconductors. The pairing amplitudes vanishes linearly at the 
first transition corresponding to the exponent Pbcs = 1 in Eq. ^80) . The Brillouin zone location 
of the gapless Dirac points in the d^2_y2 superconductor is indicated by filled circles. For the dis- 
persion Ek appropriate to the cuprates, the d^y superconductor is fully gapped, and so the second 
transition is ordinary Ising. 



low, the critical theory of the d^^-y^ to d^'^-y^ + idxy transition is not the usual (j)^ 
field theory which describes the ordinary Ising transition in three spacetime dimen- 
sions. (For the dispersion appropriate to the cuprates, the d^y superconductor is 
fully gapped, and so the d,j.2_y2 + id^y to d^y transition in Fig. 15 will be ordinary 
Ising.) 

Near the phase boundary from dr^2_y2 to dx2_y2 + idxy superconductivity it is 
clear that we can identify 

(78) 



(in the gauge where Ax2_y2 is real). We can now expand Ebcs in Eq. (75 1 for 



small (j) (with Ax2-y2 finite) and find a series with the structure |91 92 1 

Ebcs = Eq 



s4? + v\4>f 



(79) 



where s, v are coefficients and the ellipses represent regular higher order terms in 
even powers of (/>; s can have either sign, whereas v is always positive. Notice the 
non-analytic term that appears in the BCS theory - this arises from an infrared 



singularity in the integral in Eq. (75 1 over Ek at the four nodal points of the c?j.2_j^2 



superconductor, and is a preliminary indication that the transition differs from that 
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in the ordinary Ising model, and that the Dirac fermions play a central role. We can 
optimize (j) by minimizing Ebcs in Eq. (79 1 - this shows that (0) = for s > 0, 
and ((/)) ^ for s < 0. So s ~ {J2/ Ji)c — J2/ Ji where {J2/ Ji)c is the first critical 
value in Fig. [15] Near this critical point, we find 



{(/)) {Sc - S) 



(80) 



where we have allowed for the fact that fluctuation corrections will shift the critical 
point from s = to s = Sc. The present BCS theory yields the exponent Pbcs = 1; 
this differs from the usual mean-field exponent Pmf = 1/2, and this is of course 
due to the non-analytic term in Eq. (79 1. 



We have laid much of the ground work for the required field theory of the onset 



of dxy order in Section 3.2 In addition to the order parameter (p, the field theory 



should also involve the low -energy nodal fermions of the dx2_y2 superconductor, 
as described by S^p in Eq. (681. For the (p fluctuations, we write down the usual 



terms permitted near a phase transition with Ising symmetry, and similar to those in 
Eq. ([n|: 



I 



d^rdr 



(81) 



Note that, unlike Eq. (79 1, we do not have any non-analytic terms in the action: 
this is because we have not integrated out the low-energy Dirac fermions, and the 



terms in Eq. ( 8 1 1 are viewed as arising from high-energy fermions away from the 
nodal points. Finally, we need to couple the and if'i 2 excitations. Their coupling 



is already contained in the last term in Eq. (74 1: expressing this in terms of the 



fermions using Eq. ( 67 1 we obtain 



^ xy 



d^rdr 



(82) 



where dxy is a coupling constant. This coupUng also has been obtained by symmetry 
considerations, by examining invariants under the transformations of Table [T] The 
partition function of the full theory is now 



Z = j V(t)V^{D^2 exp (-5^ - 50 - S^^) , 



(83) 



where Sij, was in Eq. ( [68] l. It can now be checked that if we integrate out the W12 
fermions for a space time independent 0, we do indeed obtain a |0p term in the 
effective potential for </). 



We begin our analysis of Z in Eq. ( 83 1 by following the procedure of Section 3.2 



Assum e tha t the transition is described by a fixed point with -d^y — 0: then as in 
the theory for the transition would be the ordinary 0^ field theory 5^, 



3.2 



Section 

and the nodal fermions would again be innocent spectators. The scaling dimension 
of at such a fixed point is (1 + ?//)/2 (where r?/ is the anomalous order parameter 
exponent at the critical point of the ordinary three dimensional Ising model), while 
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that of ^^1^2 is 1 • Consequently, the scaling dimension of "dxy is (1 — ?7/)/2 > 0. This 
positive scaling dimension implies that -dr^y is relevant and the -dr^y = fixed point 
is unstable: the Dirac fermions are fully involved in the critical theory. 

Determining the correct critical behavior now requires a full renormalization 
group analysis of Z. This has been described in some detail in Ref. | [90) , and we 
will not reproduce the details here. The main result we need for our purposes is 
that couplings -dxy, u, vp/c and V/\/c all reach non-zero fixed point values which 
define a critical point in a new universality class. These fixed point values, and 
the corresponding critical exponents, can be determined in expansions in either 
(3 — d) 1 89 , 90 1 (where d is the spatial dimensionality) or |93| (where N is 
the number of fermion species). An important simplifying feature here is that the 
fixed point is actually relativistically invariant. Indeed the fixed point has the struc- 
ture of the so-called Higgs-Yukawa model which has been studied extensively in the 
particle physics literature | [94) in a different physical context: quantum Monte Carlo 
simulation of this model also exist | |95J , and provide probably the most accurate 
estimate of the exponents. 

The non-trivial fixed point has strong implications for the correlations of the 
Bogoliubov fermions. The fermion correlation function Gi = obeys 

nil \ uj + vpkxT'' + vat'' 

at low frequencies for s> Sc- Away from the critical point in the dr^2_y2 supercon- 



ductor with s > Sc, Eq ( 84 1 holds with rjj =0, and this is the BCS result, with sharp 



quasi-particle poles in the Green's function. At the critical point s = Sc Eq. (84i 
holds with the fixed point values for the velocities (which satisfy vp — va = c) 
and with the anomalous dimension 77/ 7^ - the (3 — d) expansion |89 | estimate is 



rjf Ki [3 ~ d)/14:, and the expansion estimate |93| is 77/ « l/(37r^A^), with 
N = 2. This is clearly non-BCS behavior, and the fermionic quasiparticle pole in 
the spectral function has been replaced by a branch-cut representing the continuum 
of critical excitations. The corrections to BCS extend also to correlations of the Ising 
order 0: its expectation value vanishes as Eq. (jSOjl with the Monte Carlo estimate 
/3 ~ 0.877 [95 J . The critical point correlators of 4> have the anomalous dimension 
r] w 0.754 | |95| , which is clearly different from the very small value of the exponent 
77/ at the unstable -d^y = fixed point. The value of (3 is related to 77 by the usual 
scaling law [3 ~ {1 + ri)v/2, with v w 1.00 the correlation length exponent (which 
also differs from the exponent vi of the Ising model). 



3.3.2 Nematic ordering 

We now turn, as promised, to the case of Ising-nematic ordering within the d-wave 
superconductor at a: = Xm- 

The ingredients of such an ordering are actually already present in our simple 
review of BCS theory in Section [3T| In Eq. (64 1, we introduce 2 variational pairing 
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Fig. 16 Phase diagram of Ising nematic ordering in a d-wave superconductor as a function of the 
coupling s in cS<^ . The filled circles indicate the location of the gapless fermionic excitations in the 
Brillouin zone. The two choices for s < Sc are selected by the sign of (0). 



amplitudes and Ay. Subsequently, we assumed that the minimization of the en- 
ergy led to a solution with dx2_y2 pairing symmetry with A^ = —Ay = Aj.2_y2. 
However, it is possible that upon including the full details of the microscopic inter- 
actions we are led to a minimum where the optimal solution also has a small amount 
of s-wave pairing. Then \Ax\ ^ \Ay\, and we would expect all physical properties 



to have distinct dependencies on the x and y coordinates. So, as in Eqs. (39 1, (43 i 



and ( 70 1, we can define the Ising-nematic order parameter by 



\A, 



I A, 



(85) 



The derivation of the field theory for this transition follows closely our presen- 



tation in Section 3.3.1 We allow for small Ising-nematic ordering by introducing a 
scalar field </> and writing 



A^ = A^ 



(86) 



note that I cx 4>. The evolution of the Dirac fermion spectrum under such a change is 
indicated in Fig. 16 We now develop an effective action for (f> and the Dirac fermions 
If"! 2- The result is essentially identical to that in Section 3.3.1 apart from a change 
in the structure of the Yukawa coupling. Thus we obtain a theory Sxp + 



S. 



defined by Eqs. ( 68 i and ( 8 1 1, and where Eq. ( 82 1 is now replaced by 



d rdr 



(87) 



Not surprisingly, the fermion bilinear coupling to the nematic order parameter is 
identical to that in Eq. (72 1, as is expected from the transformations of Table [T] 

The seemingly innocuous change between Eqs. ( 82 1 and ( 87 1 however has strong 
consequences. This is partly linked to the fact with 5^/0 cannot be relativistically 
invariant even after all velocities are adjusted to be equal. A weak-coupling renor- 



Quantum phase transitions of antiferromagnets and tlie cuprate superconductors 



43 



malization group analysis in powers of the coupling i9/ was performed in (3 — d) 
dimensions in Refs. 1 89 90), and led to flows to strong coupling with no accessible 
fixed point: thus no firm conclusions on the nature of the critical theory were drawn. 

This problem remained unsolved until the recent works of Refs. [88 96|. It is 
essential that the coupling -dj not be used as a perturbative expansion parameter. 
This is because it leads to strongly non-analytic changes in the structure of the (p 
propagator, which have to be included at all stages. In a model with fermion fla- 
vors, the expansion does avoid any expansion in The renormalization group 
analysis has to be carried out within the context of the 1/iV expansion, and this in- 
volves some rather technical analysis which is explained in Ref. |96 1. In the end, an 
asymptotically exact description of the vicinity of the critical point was obtained. It 
was found that the velocity ratio vp/vA diverged logarithmically with energy scale, 
leading to strongly anisotropic 'arc-like' spectra for the Dirac fermions. Associated 
singularities in the thermal conductivity have also been computed 197). 



4 Metals 

We finally turn to the transition in the metal at Xm, which anchored our discussion 
of the cuprate phase diagram in Section[T] This controls the high field transition line 
in Fig. [T] between the large Fermi surface and small Fermi pocket states. We also 
argued that this transition was a key ingredient in a theory of the strange metal. 

In addition to the order parameter ingredients we met in Sectionl2] we now have 
to also account for fermion excitations as in Section |3] In Section |3] the fermionic 
excitations had vanishing energy only at isolated nodal points in the Brillouin zone: 



see Fig. 13 In the present section we are dealing with metals, which have fermionic 
excitations with vanishing energy along an entire line in the Brillouin zone. Thus 
we can expect them to have an even stronger effect on the critical theory. This will 
indeed be the case, and we will be led to problems with a far more complex structure. 
Unlike the situation in insulators and d-wave superconductors, many basic issues 
associated with ordering transitions in two dimensional metals have not been fully 
resolved. The problem remains one of active research and is being addressed by 
many different approaches. 

As discussed in Sections |2] and 3.2 we can describe magnetic ordering by using 



either vector or spinor variables for the order parameter, and these lead to very 
different phases and critical points. For metals, the relationship between these two 
approaches, and their distinct physical properties have been described recently in 
Ref. p3) . The spinor route is more 'exotic' and leads to intermediate non-Fermi 
liquid critical phases between the small and large Fermi surface Fermi liquid phases. 
These intermediate critical phases could well be important for the experiments and 
for Fig.[T[ but we will not describe them here. We will limit our present discussion 
to the more conventional vector mode description of the SDW ordering transition. 

In recent papers 1 98 99 1 Metlitski and the author have argued that the problem of 
symmetry breaking transitions in two-dimensional metals is strongly coupled, and 
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-4. Increasing SDW order 

Fig. 17 Evolution of the Fermi surface of tlie liole doped cuprates in a conventional SDW theory 
j8| as a function of the magnitude of the SDW order \ipa\, obtained from Eq. l |90| . The right panel 
is the large Fermi surface state with no SDW order, with states contiguous to fe = occupied by 
electrons. The "hot spots" are indicated by the filled circles in the second panel from the right. 
The onset of SDW order induces the formation of electron and hole pockets (the hole pockets 
are the ones intersecting the diagonals of the Brillouin zone). With further increase of Ifa], the 
electron pockets disappear and only hole pockets remain (the converse happens in the last step for 
the electron-doped cuprates). 




proposed field theories and scaling structures for the vicinity of the critical point. 
Here we will be satisfied with a simple description of the effective action and its 
mean field theory 1 100 1: the reader is referred to the recent papers 1 98|99 1 for further 
analyses. 

As in Section|3] let us begin by a description of the non-critical fermionic sector, 
before its coupling to the order parameter fluctuations. We use the band structure 
describing the cuprates in the over-doped region, well away from the Mott insulator. 



Here the electrons Cka are described by the kinetic energy in Eq. ( 62 1, which we 
write in the following action 



(88) 



This band structure leads to the Fermi surface shown in the right-most panel of 
Fig. [17] and also later in Fig.fTS] 



4.1 SDW ordering 



As discussed in Section [3] we must now couple the fermions of Eq. (88i to the 
bosonic modes associated with the SDW ordering transition. As noted above, we 
will use the more conventional vector mode description of the SDW ordering tran- 



sition using the order parameters in Eq. ( 69 1. The analog of the coupling in Eq. ( 82 1 
now leads to the interaction term 
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Sc^ = J dTj2^^aiq)cl^K^+q,a^af!CkP+C.C.+ X ~^ y , (89) 
k.q 

where q is a small momentum associated with a long-wavelength SDW fluctuation, 
while the sum over the momentum k extends over the entire Brillouin zone. The 
complete theory for the SDW transition is now contained, in principle, in Sc + 
Sc^ + S4,, with contained in Eq. ( [7T] i. 

Let us consider the mean-field predictions of this theory in the SDW ordered 
state. For simplicity, we consider ordering at = (tTjTt), in which case ^j.^ and 



(l>ya both reduce to the Neel order ipa in Eq. ( 12 1. In the state with SDW order, we 



can take ipa = (0, 0, ip) a constant. Then Sc + iSc<j> is a bilinear in the fermions and 



can be diagonalized to yield a fermion band structure (the analog of Eq. ( 66 1) 




Ek = ± +r\ ■ (90) 



Filling the lowest energy bands of this dispersion leads to the Fermi surface structure 
|[8| shown in Fig. 17 The second panel from the right shows the Fermi surface 



obtained by translating the original Fermi surface by K, and the remaining panels 
show the consequences of mixing between the states at momentum k and k + K. 
Note that the Fermi surface has split apart into "small" electron and hole pockets, as 
discussed in Section[T] 

Let us now attempt to move beyond this simple mean field theory. As written, the 
action Sc + Sc<p + 5<f, is not conducive to a field-theoretic analysis: this is mainly 



because the sum over k in Eq. ( 89 1 extends over the entire Brillouin zone, and there 
are low-energy fermionic excitations along an entire line of k close the the Fermi 
surface. However, one simplifying feature here is that most of these low-energy 
fermions do not couple efficiently to the SDW order parameter, and their situation 



is similar to the fate of the Dirac fermions illustrated in Fig. 13 - upon scattering 
with the wavevector K, they end up at generic points in the Brillouin zone at which 
there are only high-energy fermionic states. There are now 8 special "hot spots" on 
the Fermi surface which do connect via the wavevector K to other spots directly on 



the Fermi surface: these are illustrated in Fig. 17 These "hot spots" are thus simi 



lar to the Dirac hot spots we met in Section 3.3 upon considering Ising transitions 
with a zero-momentum order parameter in a d-wave superconductor. However the 
present situation is more complex because we also have "cold lines" of zero energy 
fermionic excitations coming into the hot spots. 

A successful theory of the fermionic hot spots was reviewed in Section [33] A 
natural idea is to apply the same approach to the present situation with fermionic 
hot spots and cold lines. This leads to a problem of considerably complexity, 
which remains strongly coupled even within the context of the expansion: see 
Refs. ||99][10T) for further details. 
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Fig. 18 Phase diagram of Ising nematic ordering in a metal as a function of the coupling s in 5,^ . 
The Fermi surface for s > is the same as that in the right-most panel of Fig. [TT] but with the 
fc = point shifted from the center to the edge of the Brillouin zone. The interior regions are the 
occupied hole (or empty electron) states. The choice between the two quadrapolar distortions of 
the Fermi surface is determined by the sign of (0). 



4.2 Nematic ordering 



For completeness, we also consider the case of the Ising-nematic ordering in the 
presence of the large Fermi surface metal. Then we will have an Ising order param- 



eter represented by the real scalar field 0, which is described as before by Eq. ( 8 1 



Its coupling to the electrons can be deduced by symmetry considerations, and the 



most natural coupling (the analog of Eqs. ( 87 1 and ( 89 1) is 



The momentum dependent form factor is the simplest choice with changes sign 
under a; y, as is required by the symmetry properties of 0. Again, the sum over q 
is over small momenta, while that over k extends over the entire Brillouin zone. The 
theory for the nematic ordering transition is now described by 5c + 50 + Sc^,- The 
evolution of the Fermi surface as a function of the Ising coupling in 5^ is shown in 



Fig. 18 



Note that Eq. (91 1 does not have any large momentum transfer associated with 
K. Consequently, at any generic point on the Fermi surface, there can be scattering 
to other nearby low-energy fermionic excitations by long wavelength modes of 0. 
In other words, the entire Fermi surface is "hot". Thus we are faced with a third 
case of a "hot line" of fermions coupled to the critical order parameter mode of the 
transition. This case has been analyzed in Ref. i98|, where it is proposed that the 
critical point is actually described by an infinite number of 2+\ dimensional field 
theories, labeled by points on the Fermi surface. The reader is referred to Ref. [98 J 
for further results on this complex problem - a review of the main results appears in 
Ref. 1 1021 . 
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